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ABSTRACT 

The turbulent flow field inside a whirling annular seal was simulated by using SCISEAL, a three- 
dimensional computational fluid dynamics code. The rotor center described a circular synchronous 
whirl. A rotating frame transformation was used to make the problem quasi-steady. The flow field at an 
axial Reynolds number of 24 000 and a Taylor number of 6600 was simulated. The standard k-E model 
with wall functions and the low-Reynolds-number model were used to treat turbulence. An 
experimentally measured velocity field was used at the inlet boundary. Numerical predictions of the 
velocities and the stator wall pressures compared well with experimental data. Both turbulence 
models yielded nearly the same results. The capability of the SCISEAL code to analyze this complex 
flow field was demonstrated; the isotropic turbulence models performed adequately on this 
nonisotropic turbulence flow. 
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nominal clearance between stator and rotor, m 

rotor diameter, m 

rotor eccentricity, m 

turbulence production 

seal length, m 

static pressure, Pa 

nondimensional pressure, PL/(C A P) 

pressure drop across seal, (80.8 kPa for present case) 

rotor radius, m 

Reynolds number, 2 CpU m l\i 

position vectors, m 

Taylor number, [pw sh C/p.][2C/D] 1/2 

mean axial velocity, m/s 

Cartesian velocity components, in x, y, and z directions, m/s 

cylindrical velocity components in tangential, radial, and axial directions, m/s 

velocity vector, m/s 

normalized velocity components 

rotor surface tangential velocity, m/s 

Cartesian reference directions 

orbit radius, m 

dynamic viscosity, Pa-s 

fluid density, kg/m 3 

rotor whirl angular velocity, rad/s 

rotor spin angular velocity rad/s 



Subscripts: 


absolute 
rotor; relative 
stator 
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1. INTRODUCTION 

Turbomachinery seals are usually noncontacting and allow a leakage flow. Since the clearances are 
small, variations in the rotor position during operation can alter the fluid flow and the fluid reaction 
forces. The change in reaction forces can destabilize the rotor (e.g., a labyrinth seal) or provide 
stability (e.g., a damper seal). Evaluation of seal rotordynamic forces has long been a topic of interest. 
Black’s [1] treatment is one of the early efforts. Subsequent methods include the bulk flow models 
developed by Childs [2]. Later refinements include two-dimensional models [3] where averaged film 
properties are used. Recent advances in computational fluid dynamics (CFD) methods have prompted 
the development of three-dimensional Navier-Stokes (N-S) codes for seal flow and rotordynamic 
analyses Tam et al. [4] conducted a detailed study of the flow in seals to calculate the rotordynamic 
coefficients and to show the complexities of flows in whirling seals. A perturbation method based on 
finite-difference techniques was developed by Dietzen and Nordmann [5]. A finite-element 
perturbation model based on Navier-Stokes solutions was developed by Baskharone and Hensel [6,7] 
for seal rotordynamics. SCISEAL, a three-dimensional CFD code based on N-S equations and offering 
capabilities such as rotordynamics and different turbulence models, has been developed by Athavale 

etal. [8-10]. _ . , . . _ 

The current state-of-the-art codes rely on a variety of turbulence models to treat the turbulent flow 
fields that exist in turbomachinery seals. These models assume isotropic turbulence, and it is 
necessary to assess the accuracy of the results produced. Detailed flow-field measurements in seals for 
code validation are very few. Some studies deal with integrated quantities, such as rotor loads and 
rotordynamics (e.g., [11] and [12]). Laser Dopier anemometer measurements of flow and turbulence in 
annular and labyrinth seals were reported by Morrison et al. [13,14] for centered seals. The SCISEAL 
code was used to simulate these flows, and good correlation between the experiments and the 

numerical results was obtained [8]. , , . 

Recently, Thames [15], Morrison et al. [16,17], and Winslow [18] have reported velocity, pressure, 
and shear stress measurements in a synchronous whirling annular seal. This problem is important in 
rotor stability because the two most unstable flow-related modes in an annular seal are the rotor whirl 
at subsynchronous and synchronous whirl speeds. Such experimental data can be important t in 
validating CFD codes. The experiments were done on a seal with a large clearance-to-radius CIR 
ratio (0.015); the practical values are near 0.0015. The size of the laser probe volume dictated the 
clearance used in the seal. For practical CIR ratios, rigs with larger radius seals should be constructed 

and used. . , 

The present study was done with multiple motives. One was to assess the accuracy of the physical 
models used in the code for rotating frames. Another was to check the performance of the two most 
commonly used turbulence models: the standard k-z model with wall functions, and the low-Re model 
(both assume isotropic turbulence) in a rotating reference frame and on an nonisotropic turbulence 
problem. The rotor whirl solution is also used to calculate seal rotordynamic coefficients, and thus the 
accuracy of these results is directly related to seal dynamics. 


2. NUMERICAL METHODOLOGY 

The computations were performed by using SCISEAL, an advanced three-dimensional CFD code 
developed under NASA sponsorship for the flow and force analysis of a variety of turbomachinery 
seals [8-10] The code uses a pressure-based solution methodology to integrate the Navier-Stokes 
equations in the generalized body-fitted-coordinate (BFC) system. A finite-volume method is used to 
discretize the flow domain, and a collocated variable arrangement is used where all the velocity and 
scalar variables are stored at the center of each computational cell. Cartesian velocity components 
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are used as the primary velocity variables. The N— S equations are discretized and solved sequentially. 
A variation of the SIMPLEC method is used to provide the pressure-velocity coupling. High-order 
spatial and temporal discretization methods are available. A variety of turbulence models, including 
the standard and low-Re k-s models and a two-layer k - s model, are available in SCISEAL. The 
whirling seal problem was treated in this study by using the standard k-e model with wall functions. 
The experimental data include measurements of Reynolds stresses as well and could be used to 
validate models that predict nonisotropic turbulence. 


3. WHIRLING SEAL GRID TRANSFORMATION 

The whirling rotor generates a deforming or moving grid problem. A moving grid algorithm available 
in SCISEAL could be used, but such a time-dependent solution is expensive. When the whirl orbit is 
circular about the stator center, a rotating frame transformation can be made to make the problem 
quasi-steady. Figure 1 shows a cross section of an annular seal with a whirling rotor, with a spin speed 
CD*, a whirl speed Cl x , and an orbit radius 8. A coordinate frame Y'-Z' is attached to the rotor (Fig. 2) 
such that the Y' axis always lies along a line joining the stator center and the minimum-clearance 
point. The positions of the rotor and the coordinates Y'-Z' at two time instances are also shown. 
Obviously, in the Y'-Z' plane the flow domain does not change shape. Since the minimum- clearance 
point rotates at the whirl speed £l x , the frame X'—Y'—Z' also rotates about the X axis at Cl x . 

The rotor and stator wall velocities in the rotating frame are calculated as follows: Referring to 
Fig. 1, the absolute velocities of the rotor and stator walls can be written as 


VS.C = 0 

(1) 

= © ;c xF r +£i ;c xe 

(2) 


where r r is the position vector of a rotor wall point with respect to the rotor center. To get the veloci- 
ties in the rotating frame, one needs to subtract a velocity uq of the form 

Uq = Q. x x R (3) 


from all velocities, where R is the position vector with respect to the stator center. Thus, the stator 
and rotor wall velocities in the relative frame become 


and 

~hr=^x xR s 

(4) 

Noting that 

v r r = & x x r T + £l x x e - Q x x R r 

(5) 

we get 

R r = ? r + l 

(6) 


Vr.r =s (“*-Q*) xJ i r 

(7) 


Thus, we see that in the transformed plane the stator wall has a tangential velocity with respect to 
the stator center in the opposite direction to the whirl but that the rotor wall has a tangential velocity 
with respect to the rotor center. The apparent rotor spin is the difference between the rotor spin and 
whirl speeds. The tangential nature of the wall velocities corroborates our assertion that indeed in this 
frame the problem is quasi-steady. The momentum equations need to be modified to account for the 
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centrifugal and Coriolis accelerations resulting from the frame transformation. For a synchronous whirl 
(D x = D. x and the rotor surface velocity is zero. In effect, the whirling rotor problem has been reduced 

to a bearing problem but with the important difference that the rotor is stationary and the stator is 
spinning. This fact is crucial in understanding the flow physics and will be used in a later section to 
explain the observed results. 


4. FLOW AND BOUNDARY CONDITIONS 

Turbulent flow solutions were obtained for one of the flow conditions considered in the experimental 
data [15-16]. The experimental rig and procedures are described in [15-16]. The seal dimensions were 
L = 37.3 mm, R = 82.05 mm, and clearance C = 1.27 mm. The whirl orbit radius £ was one-half of C. 
Two different grids were used in the computations. The numbers of axial and tangential cells were 
40 and 20, respectively. In the radial direction 16 and 30 cells were used for the standard and low-Re 
models, respectively. A higher number of cells in the 0 direction would have been preferable, but the 
experimental data were available at only 20 tangential stations and hence 20 cells were used to 
expedite the interpolation of the experimental data at the inlet and exit boundaries. 

The working fluid was water with density p = 996 kg/m3 and dynamic viscosity |i = 7.83X10- 4 Pa-s. 
The fluid was taken as incompressible and with constant properties. Convective fluxes were 
discretized using the central differencing method with 20% damping (upwinding). Results presented 
here are for a rotor spin and whirl speed of 3600 rpm, which corresponds to a Taylor number of 6600. 
The mean axial velocity was 7.4 m/s, which corresponds to a Reynolds number of 24 000. 

Experimentally measured profiles of u r , uq, u x , and k were imposed at the seal inlet. A cubic spline 
interpolation procedure was used to curve fit the experimental data in the radial direction. The 
experimental data had to be modified because they were sparse and because the point nearest the 
wall was still quite far from the wall. Periodic boundary conditions were used in the 0 direction. Static 
pressure was specified at the exit boundary. Experimental distribution of static pressure in the 0 
direction was used. Transformed wall boundary velocities were imposed on the rotor and stator walls. 
As described earlier the rotor wall velocities were zero, and the stator velocities corresponded to a 
spin at a rotational speed of -Q x . 


5. RESULTS AND DISCUSSION 

Results presented here are the stator wall pressures and the contours of normalized axial, radial, and 
tangential velocities at several cross sections along the seal length. The axial u x and radial u r 
velocities were normalized with the mean axial velocity U m , and the tangential velocity uq was 
normalized with the rotor surface tangential velocity w s h- 

Figure 3 shows the distribution of the axial velocities computed by using the standard fc-e model. 
The corresponding experimental data are shown in Fig. 4. The inlet plane shows a peak velocity 
region on the suction side of the seal. This is physically consistent, since plenum flow will accelerate 
along the highest negative axial pressure gradient on the suction side. With increasing axial distance 
the peak axial velocity magnitude drops, and the peak region spreads in the larger clearance area 
(X/L= 0.49). Further downstream the peak velocity increases again, and the peak region contracts and 
moves toward the pressure side of the seal. The calculated axial velocities are in good agreement with 
the data, in terms of both the magnitude and location of the peak velocity. The numerical results show 
a thicker boundary layer in the early portion of the seal; this could be due in part to the coarseness of 
the grids, the effect of inlet boundary, and the turbulence model. The results for the corresponding 
velocities with the low-Re model are shown in Fig. 5. These results are very close to the numerical 
predictions shown in Fig. 3. 

At this point it should be emphasized again that this behavior must be explained in the transformed 
coordinates. In the absolute frame the shift of the maximum velocity region from the suction side (i.e., 
low-pressure zone) to the high-pressure side seems wrong. In addition, it is taking place in a direction 
opposite to that of the rotor spin and whirl in the absolute frame. This type of reasoning can easily 
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lead to an erroneous explanation of the results as outlined in [15-17]. Instead, one must look at the 
flow in the transformed frame, where the pressure gradients in the 0 direction are a result of the 
Couette flow generated by the moving stator wall , which “drags” the maximum velocity fluid pocket 
in a direction opposite to the rotor spin and whirl (Fig. 3). 

The normalized radial velocities for the standard k-z model are shown in Fig. 6 and the 
experimental data in Fig. 7. The inlet plane shows some effects of the vena contracta. The radial 
velocities later drop to a few percent of the mean axial velocity. The double-lobed structure in the 
latter part of the seal may be due in part to conversion from relative to absolute velocities. The results 
obtained with the low-Re model (Fig. 8) are similar to the standard k-s model results. 

The contours of the absolute tangential velocity are shown in Fig. 9. The inlet plane shows smaller 
incoming swirl that gradually increases with axial distance (spreading contours) as the rotor spins the 
fluid up. The results show good correlation with the experimental plots in Fig. 10. The low-Re results 
(Fig. 11) are similar to the results in Fig. 9. The only difference is in the contours at the seal exit, 
where the low-Re model predicts somewhat higher tangential velocities. 

The pressures on the stator wall were reported by Winslow [18], and pressure profiles from this data 
set were used as the exit boundary condition in the computations. The pressure on the outer wall as a 
function of the “time fraction” and the axial distance were measured and plots are shown in Fig. 12. 
The time fraction represents the circumferential location of a point on the stator wall when working 
with the transformed coordinates. The time fraction 0.5 corresponds to the minimum clearance; the 
pressure side (ahead of the rotor) is between 0 and 0.5 and the suction side is between 0.5 and 1.0. A 
high- pressure zone exists in the entrance region ahead of the rotor (pressure side), and a low-pressure 
zone exists behind the rotor. Along the seal length these zones reduce in intensity and get more 
evenly distributed. Near the exit of the seal the measured pressures reverse locations (i.e. the high- 
pressure zone comes over to the suction side of the rotor and the lower pressure is on the pressure 
side). In the numerical simulations the relative locations of the maximum and minimum pressures are 
maintained, although the difference between the two decreases. The experimental data show several 
“bands” along the axial directions. This feature is absent in the simulations and may be due to 
imperfections in the rotor surface. The absolute values of the maximum and minimum pressures are 
somewhat underpredicted in the simulations. Calculations show stator pressure values of -14.2 <Pc a \ c 
< 14.6 for the standard k-z model and -12.7 <P c * alc < 16.8 for the low-Re model, both of which are in 
fair agreement with the experimental values, -24.6 < P e * xpt < 19.4. 


6. SUMMARY OF RESULTS 

The turbulent flow fields in a synchronously whirling annular seal were simulated by using a three- 
dimensional computational fluid dynamics code. A transformation to a rotating frame was done to 
render the flow quasi-steady. In this frame the computed flow fields showed a fair to good agreement 
with the experimental data and a consistent picture of the flow field and physics. Although the 
turbulence models used in this study assume isotropic turbulence, the models did a good job of 
predicting flow behavior, which clearly has anisotropic turbulence, both in qualitative and quantitative 
terms. Additional flow cases with other reported flow conditions and with nonisotropic turbulence 
models are clearly needed. 
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Time fraction 



Figure 12. — Stator wall pressures as a function of axial and tangential distance. 
Time fraction 0.5 is at minimum clearance, 0 to 0.5 is pressure side, and 0.5 
to 1 .0 is suction side, (a) Experimental data, (b) Standard /c-e model results. 

(c) Low-Re model results. 
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